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Abstract 

We employ QCD sum rules to calculate the coupling constant gp^-y by study- 
ing the three point p(T7-correlation function. Our results is consistent with 
the value of this coupling constant obtained using vector meson dominance of 
the electromagnetic current and the experimental /o'^-photoproduction data. 
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I. INTRODUCTION 



The existence of the scalar-isoscalar a meson as a broad vrvr resonance has long been 
controversial. Recently, an increasing number of theoretical and experimental analyzes point 
toward the existence of this important meson Most of these analyzes find a a-pole 
position near 500-1250 MeV 0. A direct experimental evidence seems to emerge from the 
— >■ (TTT^ Stt decay channel observed by the Fermilab E791 collaboration, where a 
meson is seen as a clear dominant peak with Mo-=478 MeV and ro-=324 MeV Since a 
meson is a relevant hadronic degree of freedom, it must be incorporated into the analysis of 
hadronic processes. 

Although at sufficiently high energies and low momentum transfers electromagnetic pro- 
duction of vector mesons on nucleon targets has been explained by Pomeron exchange mod- 
els, at low energies near threshold scalar and pseudoscalar meson exchange mechanisms 
become important [Q. In particular, the presently available data on the photoproduction 
of p° meson on proton targets near threshold can be described by a simple one-meson ex- 
change model 0. In this picture, p° meson photoproduction cross section on protons is 
given mainly by a meson exchange. An important physical input for these studies is the 
coupling constant gp^-y- This coupling constant was estimated by calculating the p(T7- vertex 
assuming vector meson dominance of the electromagnetic current and then performing a 
fit to the experimental p°-photoproduction data. When this result was derived using an 
effective Lagrangian for the pcr7-vertex the value gpo-7=2.71 is obtained for this coupling 
constant Since it plays an important role in hadronic processes involving a meson, it 
is of interest to study this coupling constant from a different perspective other than vector 
meson dominance. 

In this work, we estimate the coupling constant gp^^ by employing QCD sum rules which 
provide an efficient method to study many hadronic observables, such as decay constants 
and form factors, in terms of nonperturbative contributions proportional to the quark and 
gluon condensates In section 2, we give a QCD sum rules analysis of scalar current. 
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derive a sum rule for the overlap amplitude A^- of the scalar current to the a meson state 
and estimate this amplitude since it is not available experimentally. In section 3, we derive 
a sum rule for the coupling constant gp^^ and by utilizing the experimental value of Ap, 
the calculated value of A^- and the known values of condensates we estimate this coupling 
constant. Let us here note that the nature of a meson, whether it is a qq state or it is of 
another quark, gluon structure is still a matter of debate IUJI,^. 



II. QCD SUM RULES ANALYSIS OF SCALAR CURRENT 



The QCD sum rules approach is a model independent method to study the proper- 
ties of hadrons through correlation functions of appropriately chosen currents. In the SU(2) 
flavour limit m^ = = m^, we study the scalar-isoscalar mesons by considering the scalar 
current correlation function 

n(p2) = ^ J d'^xe'P-^ < 0|r{j,(x)j,(0)}|0 > (1) 

where js = h(uu + dd) is the interpolating current for a meson. The two-loop expression for 



the scalar current correlation function n(p^) in perturbative QCD was calculated |]10[, and 
it is given by the expression 
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QCD- vacuum condensate contributions to the scalar current correlation function n(p^) were 



obtained by the operator product method and they were found in the m„=0 limit as 



n(p^ = -Q'^)cond = ^ < ^1^^ > + ^g^g2 < "^"^^ > ~27^ ^ (^s{,qqf > ■ (3) 

We note that the term < rriqqq > is not considered in the order m^ terms, as this conden- 
sate's magnitude, which is obtained by making use of Gell-Mann-Oakes-Renner relation as 
— is independent of quark mass 0. 

On the other hand, the correlation function n(p^) satisfies the standard subtracted dis- 
persion relation 



Up^rtip') = I ,(,_p2) /^(^)+n(0) (4) 

where the spectral density function is given as p{s) = (l/7r)/mn(s). We parameterize the 
spectral density as a single sharp pole n\„6{s — ml.), where the overlap amphtude is defined 
as Ao- =< 0|js|cr >, plus a smooth continuum representing the higher states. The continuum 
contribution to the spectral density is estimated in the form p = ph{s)9{s — sq) where Sq 
denote the continuum threshold and ph is given by the expression ph{s) = {1 /it) I mIlopE{s) 
with IIope{s) obtained from Eq. (2) and Eq. (3) as IIope{s) = Upertis) + Ucondis). After 
performing the Borel transformation we obtain the following QCD sum rule for the overlap 
amplitude Xa- of the scalar current 
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+2^ < -.^^ > < > < «^(^^)' > • (5) 

For the numerical evaluation of the above sum rule, we use the values < rriqqq >= (—0.82 ± 
0.10) X 10-^ GeK^ < >= (0.038 ±0.011) GeV\ < as{qqf >= -0.18 x 10^'^ GeV^ 

|l8|Jl^, and ma.=0.5 GeV. For the continuum threshold we choose So=l.l, 1.2, 1.3 GeV^ and 
we study the dependence between 0.5 GeV^ and 1.4 GeV^. The overlap amplitude as 
a function of for different values of Sq is shown in Fig. 1 from which by choosing the 
middle value = 0.9 GeV^ in its interval of variation we obtain the overlap amplitude 
as Xcr = (0.12 ± 0.01) GeV^ where we include the uncertainty due to the variation of the 
continuum threshold and the Borel parameter M^. Other sources of uncertainty are the 
errors attached to the estimated values of condensates as quoted above. If we take these 
errors into account in our analysis as well, we then obtain the value Ao- = (0.12 ±0.03) GeV"^ 
for the overlap amplitude. 



III. QCD SUM RULES FOR THE pa7- VERTEX FUNCTION 



In order to derive the QCD sum rule for the coupling constant gpa-y, we consider the 
three point correlation function 
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T,uip,p') = I d'xd^ye'^'-ye-'^-^ < 0|T{j;(0)j,^(x)j.(i/)}|0 > . (6) 

The interpolating current for p meson is = \{u'-)yU — d'j^d), and = u'yiyU + d'fiyd, 
js = ^ (uu — dd) are the electromagnetic and scalar currents, respectively. 

The theoretical part of the sum rule is obtained by calculating the perturbative contri- 
bution and the power corrections from operators of different dimensions to the three point 
correlation function T^j,. For the perturbative contribution we consider the lowest order 
bare-loop diagram. Furthermore, we consider the power corrections from the operators of 
different dimensions, < qq >, < a ■ G > and < (qq)"^ >■ We do not consider the gluon 
condensate contribution proportional to < > since it is estimated to be negligible for 
light quark systems. We perform the calculations of the power corrections in the fixed point 
gauge |jl3[. We work in the limit mg=0, and in this limit perturbative bare-loop diagram 



does not make any contribution. Moreover, in this limit only operators of dimensions d=3 
and d=5 make contributions that are proportional to < gg > and < a ■ G >, respectively. 
The relevant Feynman diagrams for the power corrections are shown in Fig. 2. 

In order to obtain the phenomenological part, we consider a double dispersion relation 
for the vertex function T^j^, and we saturate this dispersion relation by the lowest lying 
meson states in the vector and scalar channels. We obtain this way for the physical part 



Tf,u{p,P) = ^^ + ••• iV 



where the contributions from the higher states and the continuum is denoted by dots. In 
this expression, A^^ =< cr|js|0 > has been determined in Section 2. The overlap amplitude 
Xp of rho meson is defined as < 0|j^|p >= XpU^, where u^, is the polarization vector of the p 
meson. The matrix element of the electromagnetic current is given as 

< Pip)\j^Wip') >= — 9pa^K{q^){p -k Up-u-k pp) (8) 

Trip 

where q = p — p' and K{q^) is a form factor with K(0)=1. This expression defines the 
coupling constant gp^^ through the effective Lagrangian 
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C = :^gpa,dy^{d^Ap - dpA^)a (9) 



describing the p(J7- vertex 0. 

After performing double Borel transform with respect to the variables = —p^ and 
Q'^ = — p'^, and by considering the gauge-invariant structure {Pfj.qu —p - Qg/iu) we obtain the 
sum rule for the coupling constant gp^-y 

Qoa-y = — -^CM^CM'^ <UU> \ 1 771^ TTT-n o (10) 

where we use the relation < a ■ G >= < qq >. For the numerical evaluation of the 
sum rule we use the values mg = (0.8 ± 0.02) GeV^, < uu >= (-0.014 ± 0.002) GeV^ 



Tjl, and rup = 0.77 GeV, = 0.5 GeV. For the overlap amplitude \^ we use 



the value Ao- = (0.12 ± 0.03) GeV^ that we have estimated in Section 2 and for Ap we 
use its experimental value Ap = (0.107 ± 0.003) GeV^ as obtained from the measured 
leptonic width r(p° — > e+e^) of p meson |]15|. In order to analyze the dependence of gpo-7 on 
Borel parameters and M'^, we study the independent variations of and M'^ in the 
interval 0.6 GeV^ < M'^,M''^ < 1.4 GeV^ as these limits determine the allowed interval 
for the vector channel |T^. The variation of the coupling constant gpo-^ as a function of 
Borel parameters for different values of M'^ is shown in Fig. 3, examination of which 
indicates that it is quite stable with these reasonable variations of and M'^. We choose 
the middle value = 1 GeV^ for the Borel parameter in its interval of variation and 
obtain the coupling constant gpa-7 as gp^-^ = 3.2 ± 0.2 where only the error arising from 
the numerical analysis of the sum rule is considered. The other sources contributing to 
the uncertainty in the coupling constant besides those due to variations of and M'^ 
are the uncertainties in the estimated values of the vacuum condensates. If we take these 
uncertainties into account, we then obtain the coupling constant gpo-7 as gpo-7 = 3.2 ± 0.6. 

Our estimate of the coupling constant gpo-7 is in agreement with its value deduced from 
the analysis using the vector meson dominance of the electromagnetic current for the pa'j- 



vertex. On the other hand, an independent anaysis |17| utilizing the experimental value of 
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the decay rate r(p° — > 7r''"7r~7) and using the values for the a meson parameters mo-=478 
MeV and ro-=324 MeV that are obtained from the Fermilab E791 experiment gives the value 
gp(T7 = 5.92 ± 1.34 for this coupling constant. However, we like to note that in the present 
work we consider a meson in the narrow resonance limit and do not take the finite width of 
a meson into account. 



6 



ACKNOWLEDGMENTS 

We thank T. M. Aliev for helpful discussions during the course of our work. 



7 



REFERENCES 

[1] N. A. Tornqvist, Invited talk at the conference: "Possible Existence of the Light a 
Resonance and its Implications to Hadron Physics", Kyoto, Japan, 11-14 June 2000, 
KEK-Proceedings 2000-4, Dec. 2000, 224-231 and Soryushiron Kenkyu (Kyoto) 102, 
E224-E231 (2001), ( |hep-ph/0008T35| ). 

[2] N. A. Tornqvist, A. D. Polasa, invited talk at the conference "Heavy Quark at Fixed 
Target", Rio de Janeiro, Brazil, October 9-12, 2000, to appear in Frascati Physics Series, 
( |hep-ph/0011107|) . 



[3] 
[4] 
[5] 
[6] 
[7 



[9 
[10 

[11 
[12 

[13 

[14 
[15 



E. M. Aitala et al. (E791 collaboration), Phys. Rev. Lett. 86 (2001) 770. 

A. I. Titov, T. -S. H. Lee, H. Toki and O. Streltsova, Phys. Rev. C 60 (1999) 035205. 

B. Friman and M. Soyeur, Nucl. Phys. A 600 (1996) 477. 

M. A. Shifman, A. I. Vainstein, V. I. Zakharov, Nucl. Phys.B 147 (1979) 385 and 448. 

L. J. Reinders, H. R. Rubinstein, S. Yazaki, Phys. Rep. 127 (1985) 1. 

P. Colangelo, A. Khodjamirian, |hep-ph /0010175|, to be published in the Boris loffe 



Festschrift, World Scientific, Singapore, 20001. 

V. Ehas, A. H. Fariborz, Fang Shi, T. G. Steele, Nucl. Phys. A 633 (1998) 279. 

K. G. Chetyrkin, Phys. Lett. B 390 (1997) 309. 

E. Bagan, J. I. La Torre, P. Pascual, Z. Physics. C 32 (1986) 43. 

B. L. loffe, A. V. Smilga, Nucl. Phys. B 232 (1984) 109. 

A. V. Smilga, Sov. J. Nucl. Phys. 35 (1982) 271. 

V. M. Belyaev, B. L. loffe, Sov. Phys. JETP 56 (1982) 493. 

D. E. Groom et al., Eur. Phys. J. C15 (2000) 1; 



8 



[16] V. L. Eletsky, B. L. loffe, Ya. I. Kogan, Phys. Lett.B 122 (1983) 423. 
[17] A. Gokalp, O. Yilmaz, Phys. Rev. D 62 (2000) 093018. 



9 



Figure Captions: 

Figure 1: The overlap amplitude as a function of Borel parameter M^. 

Figure 2: Feynman Diagrams for the ^0-7- vertex: a- bare loop diagram, b- d=3 operator 
corrections, c- d=5 operator corrections. The dotted hues denote gluons. 

Figure 3: the couphng constant Qp^-y as a function of the Borel parameter for different 
values of M'^. 
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